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Order and disorder in irreversible decay processes 
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Dynamical disorder motivates fluctuating rate coefficients in phenomenological, mass-action rate equations. 
The reaction order in these rate equations is the fixed exponent controlling the dependence of the rate on 
the number of species. Here we clarify the relationship between these notions of (dis)order in irreversible 
decay, nA —)■ B, n = 1,2,3,..., by extending a theoretical measure of fluctuations in the rate coefficient. 
The measure, Jn — > 0, is the magnitude of the inequality between Jn, the time-integrated square of 

the rate coefficient multiplied by the time interval of interest, and £/, the square of the time-integrated rate 
coefficient. Applying the inequality to empirical models for non-exponential relaxation, we demonstrate that 
it quantifies the cumulative deviation in a rate coefficient from a constant, and so the degree of dynamical 
disorder. The equality is a bound satisfied by traditional kinetics where a single rate constant is sufficient. 
For these models, we show how increasing the reaction order can increase or decrease dynamical disorder 
and how, in either case, the inequality — £/ > 0 can indicate the ability to deduce the reaction order in 
dynamically disordered kinetics. 

PACS numbers: 82.20.Pm, 82.20.Mj, 05.40.-Fj, 05.70.Ln 


I. INTRODUCTION 

Observing the natural world, one constant seems to be 
that systems have the capacity for chemical and physical 
change, though the mechanisms for change vary. Measur¬ 
ing rates of change is a common approach for deducing 
the mechanism of kinetic processes. Complex chemical 
reactions are an example where the rate equation gov¬ 
erning molecular transformations can be evidence of a 
particular mechanism.^ For this reason, chemical kinetics 
provides a set of methods to quantify the parameters of 
the rate equation: the rate coefficient and the order(s).^ 
These empirical parameters define the dynamical equa¬ 
tions, which are the typical mass-action rate equations, 
provided the reaction system is homogeneous - as in a 
solution phase chemical reaction under standard condi¬ 
tions with uniform concentration(s) throughout. Rate 
processes with “disorder”, however, deviate from tradi¬ 
tional rate equations.^ In disordered rate processes,^ it is 
necessary to replace the rate constant in the rate equation 
with an effective time-dependent rate coefficient. The 
origin of this deviation from classical kinetics is often 
microscopic structural and energetic changes in the host 
environment, which occur on a time scale comparable to 
the chemical change.^ 

An early example of disordered or dispersed kinet¬ 
ics is the rebinding of carbon monoxide to hemepro- 
teins following photodissociation.® Experiments showed 
the complexity of this process is due to a multitude of 
energy barriers and reaction pathways. More recent ad¬ 
vances in single-molecule techniques are bringing atten¬ 
tion to the disordered decay^ of structural correlations 
in biomolecules during electron transfer reactions® and 
enzyme catalyzed reactions.Experimental techniques 
have also inspired new theoretical approaches to disor- 
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dered kinetics^®^^® in the forced unfolding®®’^^ and spon¬ 
taneous folding®®’^® of biomolecules. Disordered kinetics 
is also of interest in the charge transfer through molec¬ 
ular assemblies because of their potential use in molec¬ 
ular electronics, solar cells, and artificial photosynthe¬ 
sis.^® These diverse applications motivate the search for 
organizing principles and theoretical frameworks for ex¬ 
tracting mechanistic information that extend chemical 
kinetics. The present work develops such a framework. 

In the examples above, “disorder” can be the result 
of an intrinsically time-dependent reaction environment 
(pure dynamical disorder, homogeneous), a spatially- 
dependent reaction environment with an underlying, un¬ 
known distribution of rate coefficients (static disorder, 
heterogeneous), or both. In any case, the overall kinet¬ 
ics is describable by a single time-dependent rate coeffi¬ 
cient, with underlying rate coefficients that can fluctuate 
in time or space. The reaction order of the overall kinet¬ 
ics, however, is a constant; as in traditional kinetics, the 
reaction “order” in the rate equation is the power of the 
concentration dependence of the rate. While the reaction 
order can be evidence of a particular mechanism in tradi¬ 
tional kinetics, it is not clear whether the reaction order 
might provide mechanistic clues in kinetics with dynami¬ 
cal disorder. As a step in this direction, our interest here 
is in the temporal variation of an effective, overall rate 
coefficient for disordered kinetic processes, and how this 
disorder relates to reaction order. 

Relaxation or irreversible decay phenomena have been 
the primary focus of theoretical and computational works 
in disordered kinetics.® There has been a particular the¬ 
oretical emphasis on first-order decay"*; however, un¬ 
der some circumstances, disordered processes can have 
higher-order rate equations; for example, this was shown 
for a case of dispersive, early stage protein folding kinet¬ 
ics.Complicating the mechanistic insights that can 
come from rates, is that experimental data may agree 
with both first- and second-order rate equations when 
the rate coefficient is time dependent. A further com- 
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plication is that multiple microscopic models can lead to 
the same macroscopic relaxation behavior; an example 
for first-order decay are the mechanisms consistent with 
the stretched-(non)-exponential decay.^^ These findings 
stimulate our investigation into how meaningful the tra¬ 
ditional notion of reaction-order is for phenomenological, 
non-exponential kinetic data. 

Here we consider the dynamic disorder in the macro¬ 
scopic kinetics of the irreversible processes 

nA^B for n= 1,2,3,.... (1) 

We assume the concentration, Cn{t), of the species A 
obey the (non-)linear differential rate equations 

= -knit) [Cnit)r (2) 

with time-dependent, effective rate coefficients knit). 
That is, the rates are n*^-order in A. [The subscript 
n will distinguish variables and parameters of rate pro¬ 
cesses of different reaction order.] Since these rate equa¬ 
tions are solvable for any n, we will consider the entire 
class of reactions. However, the highest known order of 
a chemical reaction is third order in traditional kinetics, 
kit) —)■ In disordered kinetics, second-order seems to 
be the highest known order in agreement with experi¬ 
mental results. For example, power law decays describ¬ 
ing recombination reactions and peptide folding satisfy 
second-order rate equations,though stretched expo¬ 
nential functions also fit the available data.^^’^^ General¬ 
ized functions, encompassing both stretched exponential 
and power law decay are also known. 


II. THEORETICAL BACKGROUND 


Static and dynamic disorder lead to an observed rate 
coefficient, knit), that depends on time in macroscopic 
rate equations. The time dependence can stem from fluc¬ 
tuations in the rate coefficient, fluctuations due to the 
heterogeneity of the system and/or the reaction environ¬ 
ments. Measuring the dispersion of the rate coefficients 
is therefore a natural approach to quantify the involve¬ 
ment of the environment in the mechanism or the struc¬ 
tural disorder of the medium.^® The degree of variation 
in the observed rate coefficient also quantifies the hdelity 
of a rate coefficient or rate equation, assuming rate co¬ 
efficients that vary less over the observational time scale 
are more desirable. 

A cumulative measure of disorder that we extend here, 
originally in Reference 29, is the magnitude of the in¬ 
equality 


JniAt) > CniAt)^ 


between the statistical length (squared) 


AniAt)^ 



1 


k„it)dt 


(3) 

(4) 


and the divergence 

JniAt) = At knit)‘^dt (5) 

over a time interval At = tf — ti- From these definitions, 
Ln and Jn are functions of a possibly time-dependent 
rate coefficient. The statistical length, we use is the 
cumulative time-dependent rate coefficient over a period 
of time At, and the divergence, Jn, is the cumulative 
square of the rate coefficient, multiplied by the time in¬ 
terval. These definitions are independent of the particu¬ 
lar system and could be applied to any disordered kinetic 
process with a time-dependent rate coefficient. The dif¬ 
ference, Jn — C^, quantifies both the variation of knit) 
in time and the fraction of the observation time where 
those changes occur. As a measure of disorder, the in¬ 
equality is complementary to other information theoretic 
measures of dynamic disorder.®®’®^ 

The representation of Jn and Cn we dehne above is 
just one from our earlier work on first-order (n = 1) rate 
processes,^® which established a connection between the 
rate coefficient and a modified Fisher information.®® For 
first-order rate processes, the inequality above has two 
important features: (1) Ji—C\ measures the variation of 
the rate coefficient in statically or dynamically disordered 
decay kinetics over a time interval of interest, and (2) 
the lower bound, J\ = C\, holds only when the effective 
rate coefficient is constant. These features are sensitive 
to the definition of the time-dependent rate coefficient, 
knit), but we will show that defining it appropriately, the 
inequality applies to irreversible decay kinetics with non¬ 
linear rate laws (i.e., kinetics with total reaction-order 
greater than unity). We illustrate this framework with 
proof-of-principle models for irreversible decay phenom¬ 
ena of total reaction-order greater than one. 

In macroscopic chemical kinetics, experimental data 
are typically a concentration profile corresponding to the 
integrated rate equation. For a population of species irre¬ 
versibly decaying over time, we define the survival func¬ 
tion for n*^*'-order decay. 

For truly irreversible decay, this is the probability the 
initial population survives up to a time t. In traditional 
first-order kinetics the survival probability, is char¬ 

acterized by a rate constant wi. Linear htting methods of 
survival data are a standard approach to find an empiri¬ 
cal rate equation. For example, for a first-order process, a 
plot of In Si it) versus time t is linear with a slope of — wi 
only if the traditional first-order rate equation is valid. If 
the graph is non-linear, it is natural to search for a higher- 
order rate equation that gives a linear graph.® Another 
approach is to introduce a time-dependent rate coefficient 
ki it) describing the non-exponential decay, when one sus¬ 
pects members of the population decay in different struc¬ 
tural or energetic environments or in a local environment 
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that fluctuates in time. This approach should extend to 
higher reaction orders, but seems to have received less 
attention than the case of first-order kinetics. 

We define the effective rate coefficient, through 

an appropriate time derivative of the survival function, 
depending on the reaction-order, n, 

{ —-f- In^ift) if n = 1 

“^n-Idt ^ ■ 

In any case, A:„(t) has dimensions of inverse time. The 
survival function is dimensionless, which prevents ki{t) 
from having a logarithmic argument with units. While 
the survival function is not necessary on the grounds of 
units for n > 2, we use it to keep the interpretation of 
S{t) as a survival probability and to maintain consistency 
with our previous work.^® 

With these we can express the statistical length, 
Cn, and the divergence, J^n, for any reaction order. It is 
then straightforward to show these kn{t), give the bound 
J^n = if kn{t) is independent of time. As an example, 
consider an n‘^-order reaction where n > 2. The tradi¬ 
tional integrated rate equation is the time dependence of 
the concentration, C'n(<), 

^ Cn(o)"-i ^ 

When n ^ 1, the rate constant, has dimensions of 
[time]“^ • [concentration]”*^From the general defi¬ 
nition in Equation 6, the survival function is 


Snit) = 


1 


1-b (n - l)w„C'„(0)" If 


(9) 


and the definition of the effective rate coefficient is 
kn{t) —This time-independent rate co¬ 
efficient gives for both the statistical 

length squared and the divergence. That is, the bound 
J'n = holds when there is no static or dynamic dis¬ 
order, and a single rate coefficient, a;„, is sufficient to 
characterize irreversible decay. When the kinetics are 
statically or dynamically disordered, the rate coefficient 
governing irreversible decay is not constant and one must 
work with the above definitions of knit). 


exponential decay is “stretched” to characterize non¬ 
exponential decays with a two-parameter survival func¬ 
tion exp[—(wit)^]. The parameter oji is a characteristic 
rate constant or inverse time scale and the parameter /3 
is a measure of the degree of stretching; both parame¬ 
ters depend on the system and can depend on external 
variables such as temperature and pressure.® This proce¬ 
dure for modifying traditional kinetics introduces a time- 
dependent rate coefficient, ki^pit) = l3{oJit)^/t, into the 
empiricaf modef. While different mechanisms can lead 
to stretched exponential decay,^® and any model with 
a time-dependent rate coefficient can be subject to our 
analysis, this type of decay is a good example for our 
purposes, given its generality.^"*" 

We introduce disorder and a time-dependent rate coef¬ 
ficient into “higher-order” kinetics through the n*®-order 
survival function 


Sn,0it) = 


1 -b (n - 1) [a;„C'„(0)"' 


( 10 ) 


Again there is a characteristic rate constant a;„ and dis¬ 
order parameter 0 < /3 < 1. We chose = 1 and a 
fixed observational time scale to simplify the presenta¬ 
tion of our results. Plonka and coworkers interpret the 





III. RESULTS AND DISCUSSION 
A. Kinetic model with disorder 

To study the effect of dynamical disorder on kinet¬ 
ics with reaction orders greater than one, we adapt the 
Kohlrausch-Williams-Watts (KWW) model.The 
procedure is similar to that of first-order kinetics where 


FIG. 1. When the disorder parameter, /?, is (a) 0.25, (b) 0.50, 
or (c) 0.75 the function t/S"~^^{t) is non-linear in time. This 
non-linearity motivates the use of an effective time-dependent 
rate coefficient for the overall decay when the rate equation is 
known to be n*®-order in the decaying population. Traditional 
n*’'-order kinetics corresponds to the special case of (d) /3 = 1, 
where the rate coefficient is time independent and comes from 
the slope of these plots, uj„Cn{0)"~^■ Here CniO) = 2 and 
= 1 in arbitrary units. The effective, time-dependent rate 
coefficient is proportional to the slope at a given time. 
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FIG. 2. Time-dependence of the effective rate coefficient 
kn,p{t) for varying reaction orders n = 1, 2, 3, and 4 when 
P is (a) 0.25, (b) 0.50, (c) 0.75, and (d) 1.00. The kn,p{t) of 
different order n become the same function of time as /I —l 0 
and time independent as /? —>■ 1. In all cases the initial con¬ 
centration is C'n(O) = 2.0 and the rate constant is a;„ = 1 in 
arbitrary units. 


stretching parameter for n = 2, as the result of a time- 
dependent energy barrier height, and as the result of a 
superposition of second-order decays.^® For n > 2, the 
characteristic kinetic plot - a graph of 1/5'"^^ versus t 
- is shown in Figure 1. As these graphs illustrate, de¬ 
creasing the parameter /3 stretches the survival plot for 
n*''-order kinetics, just as it does for first-order processes. 
Another similarity with first-order kinetics is that the 
limit /3 = 1 corresponds to traditional kinetics, the ab¬ 
sence of disorder, and when the effective rate coefficient, 
kn{t) —> WnCn(0)"“^, corresponds to the constant slope 
of the graph. 

At and away from the limit /3 = 1, the survival func¬ 
tions satisfy rate equations 

= -k^At)SnAtr ( 11 ) 

with the same form as Equation 2. The effective rate 
coefficients in these cases are 

kn,l3{t) = P [uJnCn{0T~^t]^ /t 

= ( 12 ) 

This definition of the time-dependent rate coefficients dif¬ 
fers from that of other authors^^ to satisfy the rate equa¬ 
tions above. Keep in mind, we define the effective rate 
coefficient in terms of the dimensionless survival function 
Sn,i 3 {t). If we instead define the effective rate coefficient 


in terms of the concentration, Cn^pit), then the effective 
rate coefficient becomes a constant as /3 —>■ 1, with di¬ 
mensions of [concentration]• [time]“^. We will use 
the former definition of the effective rate coefficient be¬ 
cause the subsequent results will depend on the initial 
concentration of the decaying species, which is consis¬ 
tent with other results in traditional kinetics (e.g., the 
half-life depends on the initial concentration in kinetics 
with n > 2).^ 


B. Inequality for disorder in irreversible kinetics 


Having laid out the details and main aspects of the 
model, we turn to the application of the theory. For this 
class of disordered kinetics, it can be difficult to distin¬ 
guish the survival curves for different n when /? is small. 
For example, experimental P values around 0.1 are nec¬ 
essary to fit peptide folding data to a stretched expo¬ 
nential.In this case, the data are also well fit by 
an asymptotic power law decay (the integrated form of 
a second-order rate equation). Our results support this 
fact, and show it is also true of irreversible decay for 
higher n. Figure 2, shows the time dependence of kn,p{t) 
for u = 1, 2, 3, and 4 when P is (a) 0.25, (b) 0.50, (c) 
0.75, and (d) 1.00. These data are for an initial concen¬ 
tration of (^^(O) = 2.0 and = 1.0 in arbitrary units. 
As /3 —)■ 0 the effective rate coefficient —)■ kpA, 

and the corresponding survival plots of all reaction or¬ 
der converge, Sn,p{t) —)■ Sp{t). That is, it can become 
difficult to distinguish the order of a reaction when P 
is small, which can be taken to mean decay events are 
highly cooperative or correlated. 

When a P approaching zero is necessary to fit the avail¬ 
able survival data, processes with an order n > 2 may not 
be distinguishable from first-order processes. That is, 
for this widespread type of kinetics, higher-order kinetics 
may be an equally valid description of experimental data 
when P is small, unless the modeller has more informa¬ 
tion about the underlying mechanism to suggest a par¬ 
ticular reaction order. The greater the value of /? is, than 
say 0.1, the easier it should be to discern rate equations 
with n > 2 within the experimental uncertainty. Exactly 
how small P must be for the reaction order to lose mean¬ 
ing will depend on the value of a;„, which we take to be 
unity, and the observational time scale. When the rate 
equations of different orders are distinguishable, we will 
show there is an interesting interplay between order, n, 
and the disorder parameter, /3, that we can disentangle 
with the measure of dynamical disorder, — Cn,/ 3 - To 
do this, we will need the length and divergence. 

From the effective rate coefficient for n > 2, the statis¬ 
tical length is 


= [^nCr,i0r-H]^ \l' 


(13) 
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FIG. 3. The magnitude of the inequality, on a 

logarithmic scale over the range of the disorder parameter, /3, 
when the initial concentration C„(0) is (a) 0.5 and (b) 2.0. 
Here ujn = 1, ti = 1, and f/ = 10 in arbitrary, but consistent, 
units. Increasing the reaction order, n, can increase or de¬ 
crease the dynamical disorder, as measured by the inequality, 
depending on the initial concentration. The inequality cor¬ 
rectly reflects the absence of dynamical disorder in the limits 
/3 0 and /3 = 1. 


and the divergence is 


Jn„s(Af) 

At 


( ^ 2 / 3-1 
X J 2/3-1 


if/?^i 

if/3=i. 


(14) 


The effective, time-dependent rate coefficient depends on 
the initial concentration, which carries through to J^n,p 
and Cn,p- Comparing the inequality for first-order kinet- 
ics^® to that of higher-order kinetics, we see the inequal¬ 
ities are related: 


- Clp = C„(0)2/5("-i) - Clp] . (15) 

The initial concentration, and its dependence upon n and 
/3, distinguishes the inequality of higher-order kinetics 
from first-order kinetics. We will examine this relation 
as a means of determining the reaction order, n, in dis¬ 
ordered kinetics after exploring the effect of n and /? on 
the inequality. 


C. Disorder and the parameters f3 and n 

From these results we can show Sn,p ~ & measures 

dynamical disorder for irreversible decay kinetics with 
n > 1 . A minimum criterion a reasonable measure of 
disorder must have is that it be zero in the absence of 
disorder. For the present model, this means the measure 
must then be zero when /3 = 1 and tend to zero as —>■ 0 , 
since the rate coefficient kn^p{t) is time independent in 
these limits. Indeed, we find the inequality satisfies this 
requirement, and the lower bound Sn,p ~ ^ = 0 holds 

when /3 = 0, 1 irrespective of the reaction-order, n. It is 


important to note that this is true when kn{t) is defined 
as in Equation 12. 

Away from the limits ,0 = 0 and 1, we find the statisti¬ 
cal length and divergence are not equal, Jn,p — ^ > 0 

(Figure 3). Interestingly, for all n, their difference ini¬ 
tially increases with decreasing /3 from one to zero, since 
this stretches the survival plot and increases the variation 
in the effective rate coefficient over the observational time 
window. However, as /? decreases further, there is even¬ 
tually a point at which the rate coefficient is effectively 
constant over the majority of the observational time win¬ 
dow. This leads to an eventual decrease in disorder when 
stretching n‘*'-order kinetics due to the separation of de¬ 
cay into “fast” and “slow” components. The inequal¬ 
ity accurately reflects this effect and, consequently, has a 
maximum in its dependence upon /3, tending to the lower 
bound Jn,p — ^ = 0 as /3 —>■ 0. Another interpreta¬ 

tion of the inequality is as a measure of the deviations of 
the 7T,**'-order characteristic kinetic plot from linear dur¬ 
ing observational time window, and the portion of the 
time window when those deviations occur. The maxi¬ 
mum in the inequality results from the competition of 
these factors. These findings are evidence that the in¬ 
equality measures the amount of dynamical disorder in 
irreversible decay processes of any reaction-order. 

Figure 3 shows the dependence of the degree of dynam¬ 
ical disorder on the disorder parameter, f3. The inequal¬ 
ity reveals a main feature of dynamical disorder: the dis¬ 
order parameter, /3 does not linearly increase the amount 
of dynamical disorder in irreversible relaxations describ- 
able by stretched exponential and power law (n > 2 , 
Equation 10) decay. Further, the amount of dynamical 
disorder measured by the inequality does not uniquely 
determine the value of f3. Increasing (for /3 values greater 
than that corresponding to the maximum in ^) 

or decreasing /3 (for (5 values less than that correspond¬ 
ing to the maximum in Jn,p — ^) can improve the 

approximation of the decay over the observational time 
scale with a single rate coefficient. We take this as fur¬ 
ther evidence that the inequality is a useful measure of 
dynamical disorder. 

Comparing panel (a) and (b) in Figure 3, we see 
the initial concentration dictates whether increasing the 
reaction-order increases or decreases the inequality, and 
so, the dynamical disorder for a given /3, a;„, and obser¬ 
vational time scale. Panel (a) shows that for an initial 
concentration C'„(0) < I, Jn,p — p decreases with re¬ 
action order, n, at a given /3 7 ^ 0 or 1. This suggests 
that for a given /3, higher-order the kinetics will have 
an effective rate coefficient kn,p{t) that varies less over 
the observational time scale: there is less dynamical dis¬ 
order with increasing reaction order. In contrast, panel 
(b) shows the converse is true when the initial concen¬ 
tration C'„(0) > 1: Jn,^ ~ ^ increases with reaction 

order, n, suggesting kn,p{t) varies more over the obser¬ 
vational time scale the higher-order the kinetics. Over¬ 
all, higher reaction-order processes can have more or less 
disorder for a given (/3, a;„. At) depending on the initial 
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C„(0) = 0.5 C„(0) = 2.0 



FIG. 4. The magnitude of the inequality, on 

a logarithmic scale over the range of the reaction-order, n, 
when the initial concentration C„(0) is (a) 0.5 and (b) 2.0. 
Here = 1, = 1, and tf = 10 in arbitrary, but consistent, 

units. The disorder parameter, /3, controls how strongly the 
dynamical disorder, as measured by the inequality, depends 
on the reaction-order, n, regardless of the initial concentra¬ 
tion. When P = 1, Jn,p = for all n. 


concentration. The reaction-order is not the opposite of 
dynamical disorder. 

Among the growing number of disordered systems,^’® a 
couple of interesting applications of our approach would 
be peptide folding^^ and low-temperature recombination 
reactions, such as the recombination of carbon monoxide 
with myoglobin,® which have small stretching parame¬ 
ters (/3 « 0.1). Inspecting Figure 3 around P = 0.1 we 
find there are n-dependent differences in the ^ 

curves. At /? = 0.1 the inequality is small and only 
weakly dependent on the reaction order. This is quanti¬ 
tative evidence of why experimental peptide folding data 
that is well fit by a stretched exponential (n = 1 ) with a 
P value around 0.1 is also fit by an asymptotic power law 
decay (u = 2). The near equivalence of j 7 i ./3 — p and 
p for reactions of different n suggests the de¬ 
gree of disorder is approximately the same within exper¬ 
imental uncertainties. This result also suggests higher- 
reaction-orders may fit this data. 

Using this model for irreversible decay, we can also con¬ 
sider the dependence of the inequality on the reaction- 
order parameter, n. Figure 4 shows that /3 deter¬ 
mines how strongly the inequality (dynamical disorder) 
depends on the reaction order. We see increasing p 
strengthens the dependence of the dynamical disorder 
on the reaction-order - whether the dynamical disorder 
increases or decreases with n. In particular, as /I —)■ 0 
the degree of dynamical disorder, as measured by the in¬ 
equality, becomes more uniform over the range of n; the 
dynamical disorder becomes independent of the reaction 
order. Furthermore, this finding is independent of the 
initial concentration, as the two panels show. We also 
saw that as /3 —> 0 , the n-dependent variations in the 
effective rate coefficient vanish, kn^p{t) —>■ 0 in Figure 2, 
which is why Jn,p — C?n ^ ^ 1 ®° becomes independent of n. 
Since the parameter P determines the rate of change of 
Jn,i3 — 13 with n, for small p, the lack of variation of 




FIG. 5. (a) A log-log plot of the inequality, Jn,p — as 

a function of the initial concentration, C'n(O), for n = 2 — 4 
and P = 0.1, 0.5. (b) The inequality, J 2 ^p — C^ p, for second- 
order kinetics on a log scale versus the disorder parameter, 
P, for initial concentrations, C'n(O), spanning four orders of 
magnitude and a fixed observational time scale {ti — 1, tf = 
10 ). 


the inequality with n suggests a fundamental difficulty in 
determining the reaction order. These results are another 
reflection of the difficulty in determining the reaction or¬ 
der in non-exponential kinetics when P is small. 


D. Disorder and the initial concentration, Cn(0) 

In traditional macroscopic kinetics, initial rate meth¬ 
ods are a standard means of determining the order(s) 
in the rate equation. Initial rate methods involve vary¬ 
ing the initial concentration of a particular reactant, and 
measuring the initial rate. Through a series of such 
measurements, an experimenter can find the order of a 
hypothesized rate equation through algebraic manipula¬ 
tion, and gain clues about the mechanism of the reac¬ 
tion. We take a similar approach for the present model 
of disordered kinetics. Figure 5(a) shows a log-log plot 
of J'n,/3 — p as a function of the initial concentration, 
C„(0), for n = 2 — 4 and P = 0.1, 0.5. The inequality is 
over the observational time window, ti = 1 and t/ = 10. 

Most immediate from these data is that the magni¬ 
tude of the inequality, and the dynamical disorder it rep¬ 
resents, directly correlates with C'„(0) for all P and n. 
At a given /3, Jn,^ ~ p increase or decrease with 
reaction-order, n, depending on the initial concentration: 
the inequality increases with n for C'„(0) > 1 and de¬ 
creases for Cn{0) < 1, as shown in Figures 3 and 4. Also 
of note is that while J'n^p ~ p increases more rapidly 
for n = 4 than n = 2, the difference can be offset by a 
small P, e.g., 0.1. This detail reflects the more general 
finding that the “strength” of the increase or decrease of 
the inequality depends on P; decreasing the value of P 
suppresses the concentration dependence of the dynam¬ 
ical disorder (inequality) for all reaction orders, and the 
differences in the inequality for different n. Thus, even 
the initial concentration dependence may not be suffi¬ 
cient to distinguish the disordered kinetics of reactions 
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with different reaction order. 

Figure 5(b), shows how the inequality for second-order 
kinetics, f72,/3 — p, varies with the disorder parameter, 
/3, for initial concentrations, a;„C'„(0), spanning four or¬ 
ders of magnitude and a fixed observational time scale. 
The maximum in the inequality at /3niax shifts to larger 
/3 values with increasing (72(0). Similar behavior re¬ 
sults from increasing uj 2 , which we have set to one. For 
/3 < /3max5 P and n have the opposite effect on the dynam¬ 
ical disorder: increasing /3 and decreasing n will increase 
the dynamical disorder. But if /3 > Pmax, P and n have 
the same effect on the dynamical disorder: decreasing 
P and decreasing n will increase the dynamical disor¬ 
der. We conclude that n and P are not strictly opposed 
parameters of order and disorder, and the inequality, a 
measure of dynamical disorder is necessary to disentangle 
these subtleties of disordered kinetics. 

Measuring the initial concentration dependence of the 
inequality is a straightforward procedure to try to de¬ 
termine the reaction order, and a clue to the underlying 
mechanism, even if the kinetics is disordered. However, 
the inequality reveals that in disordered kinetics there is 
a fundamental difficulty in identifying a unique assign¬ 
ment of the reaction order from survival data. The dif¬ 
ficulty will surely increase when there is uncertainty in 
the survival function, as will be the case in experiments 
and simulations. Under these circumstances, the fluc¬ 
tuations in the survival function will also likely impact 
the inequality, though here the concentrations and sur¬ 
vival functions are non-fluctuating quantities. Even un¬ 
der ideal conditions, as we have shown, the reaction order 
is less meaningful when a small disorder parameter, P, is 
necessary. 


IV. CONCLUSIONS 

In summary, we have considered the inequality, 
Jn > iPii, between two functions of the effective time- 
dependent rate coefficient, the statistical length and di¬ 
vergence. We have extended the scope of the inequality 
to irreversible decay processes with any reaction-order, 
by appropriately defining kn(t). A redefinition of fc„(t) 
preserves the advantageous features of the inequality. 
First, the inequality quantifies both the variation of fc„(t) 
in time and the fraction of the observation time where 
those changes occur - the amount of dynamical disorder. 
Second, the lower bound of the inequality, is 

a condition for traditional kinetics to hold, and for rate 
coefficient to be unique. Leveraging this measure, we 
have found for stretched exponential and power law de¬ 
cay that the dynamical disorder can increase or decrease 
with the reaction order. Higher reaction order does not 
imply lower dynamical disorder, and vice-versa. Instead, 
the parameter P determines the ability to distinguish ki¬ 
netics of different reaction order in the types of decay 
considered. As in traditional kinetics, the initial concen¬ 
tration may be used to distinguish reactions of different 


overall order. Here, for dynamically disordered decay, the 
initial concentration dependence of the inequality proved 
to be a useful way to determine the reaction order (when 
such a parameter is meaningful), establish a rate equa¬ 
tion, and gain clues about the reaction mechanism. From 
our results, it seems justifiable to consider higher-order, 
n > 1, processes for experimental measurements of decay, 
even when there is disorder. 
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